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We study the leading long-distance attractive force between two holes in a plate arising from 
a scalar field with Dirichlet boundary conditions on the plate. We use a formalism in which the 
interaction is governed by a non-local field theory which lives on the two holes. The interaction 
energy is proportional to QxQijr 1 at large separation r, where Q\ and Q2 are certain charges 
associated with the holes. We compute these charges for round and rectangular holes. We show 
that the 1/r 7 behavior is universal for separations large compared to the linear dimensions of the 
holes, irrespective of the spin or interactions of the bosonic field. We also study the interaction 
between two long thin slits, for which the energy falls off as 1/r 6 . 



I. INTRODUCTION 



The original Casimir effect described the interaction 
between two parallel conducting plates due to vacuum 
fluctuations of the electromagnetic field pQ. Since the 
pioneering work of Casimir many variants of this effect 
have been studied. For a recent review see [2]. 



In the present paper we consider a single infinite plate 
with two holes in it. We consider a scalar field with 
Dirichlet boundary conditions on the plate, and study the 
way in which the holes modify the ground state energy of 
the field. We do this using a formalism developed in [3J, 
where one first integrates out the scalar field in the bulk 
to obtain a non-local field theory that lives in the holes. 
This description makes it easy to study the leading long- 
distance interaction between holes. We find that the force 
between the holes is attractive and that the interaction 
energy scales as 1/r 7 , where r is the distance between 
the holes. The interaction energy is also proportional to 
the product Q1Q2 of certain charges associated with the 
holes. These charges have units of (length) 3 and depend 
on the geometry of the holes. 



An outline of this paper is as follows. In section |TT] 
we set up the general formalism and extract the distance 
dependence. We then compute the charge for round and 
rectangular holes. In section |HI| we study the interac- 
tion between two long thin parallel slits, using a mix of 



numerical and perturbative techniques. Section IV gives 



II. INTERACTION BETWEEN HOLES 

Consider a scalar field </> in the presence of an infinite 
plate with two holes in it, as shown in Fig. [l] We first 
work in three Euclidean dimensions with a scalar field of 
mass /i. We introduce a pair of coordinates x along the 
plate and a single transverse coordinate y. We may take 
the field <j> to be confined to a large cubical box with the 
plate of interest partitioning this box into a left region 
and a right region. The location of the plate is at y = 0. 
We impose a Dirichlet condition on the plate. 

The approach of [3] begins by integrating out the scalar 
field in the bulk to write a non-local effective action for 
the field in the holes. One way to think about this effec- 
tive action is to regard y as Euclidean time. The ground 
state wavefunctional for the field takes the standard form 

HE] 



^o[<P\ = const, exp — 



J d 2 x l -4>V^~ 2 



+ fi 2 <P 



(l) 



Here 4>{x) is the value of the field on a slice of fixed y 
and V 2 is the Laplacian on K 2 . The three-dimensional 
partition function of the field in the presence of the plate 
at y — is then 



z 3d = ^$;^]$oW 



(2) 



the argument for the universality of the 1/r 7 interaction. 
We conclude with some remarks in section fVl 
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FIG. 1: An infinite plate with two holes separated by a dis- 
tance r. The field vanishes at the location of the plate. In 
the holes we denote the fluctuating value of the field by <f>o- 



where <j>q is the value of the field at y = 0. The boundary 
conditions imposed by the plate are taken into account 
by requiring that 4> vanishes outside the holes. This 
formalism was used in [3] to obtain the diffractive edge 
effects for the Casimir interaction between a plate and 
half-plate parallel to it, as well as for a plate and another 
plate perpendicular to it at a finite separation. This gave 
results in agreement with numerical calculations [5] for 
both cases and with a special calculation for the case of 
perpendicular plates [7]. 

For working out the functional integral in pi explicitly 
for the present case of two holes on the plate, it is useful 
to introduce projection operators Pi , P 2 on the two holes. 
These are operators on L 2 (R 2 ), defined by 



Pif(x) = 



m 



if x G hole i 
otherwise 



The Gaussian functional integral in pi) then gives 



lo g ^Trlog(g;g : 



(3) 



(4) 



where Oij = Pi\J—V 2 + H 2 Pj- To extend this result to 



four dimensions we think of the mass /i as arising from 
Kaluza-Klein momentum around an additional periodic 
coordinate with period /3 — > oo. This means the four- 
dimensional partition function is 



log Z 4d = fi 



2n~ 



(-\ogZ 3d ) 



(5) 



The interaction between holes is clearly due to the off- 
diagonal entries Oi 2 , ©21 • If the holes are widely sepa- 
rated these off-diagonal entries will be small and we can 
treat them pcrturbatively. The zeroth order term in this 
perturbation series is independent of the separation dis- 
tance, and the first order term vanishes identically. So 
the leading contribution to the interaction energy is 



E„ 



1 f°° 
— J dfiTr [(0 11 )' 1 O l2 (0 22 )- 1 21 \ 



(6) 



In this equation, for the inverses, one first constructs the 
diagonal projected operators On, then inverts them in 
the subspace corresponding to hole i. Writing the trace 
out in position space we get 



E-„ 



1 

27T 



1 t 

dp I d 2 x 1 d 2 x' 1 d 2 x 2 d 2 x' 2 (xi| (On)" 1 |x' x ) (x' 1 |Oi 2 |x 2 ) (x 2 | (0 22 y l |x 2 ) (x 2 |G 21 |xi) 



(7) 



where xi, x' x correspond to points in the first hole and 
x 2 , x 2 to points in the second. 

For widely-separated holes there are further simplifica- 
tions we can make. Large separation necessarily means 
the separation distance r is large compared to the size of 
the holes, so we can pull the off-diagonal matrix elements 
out of the integrals. Indeed we have (the projection op- 
erators are immaterial for this) 



(xilOulxa) 



1 



fi 2 \r) 

00 ds 
S3/2 \' 



S (-V 2 +,u 2 



I''') 



27T7- 3 



(1 + fir)e 



(8) 



the explicit fi dependence gives the leading long-distance 
interaction between small holes, 



E„ 



5QiQ 2 



327r 3 r 7 
where the charge associated with each hole is 



(10) 



Qi = / d 2 xd 2 x' (x| (Cu)- 1 |x') • (11) 

./hole i 

Note that Qi is simply the matrix element of 
(PiV— V 2 Pi)^ 1 between functions that are constant 
(equal to one) in the hole. 



where we introduced an integral representation for the 
square root and used the heat kernel 



|e sV \r) 



1 

4ns 



a -r 2 /As 



(9) 



In 



princip 



le the 



remaining 



matrix elements 



(xiKOnrVi), (x 2 |(0 22 ) _1 |x 2 ) depend on M . 
But for small holes we can neglect this dependence. 
(Retaining it would give corrections down by powers of 
(size of hole)/r.) So using ^ in Q and integrating over 



Round holes 



The charge Q defined in (11 1 depends on the geometry 



of the hole in question but is independent of the separa- 
tion distance. Here we study it for a round hole of radius 
R. 

The approach we will use is quite si mple: we adopt a 
lattice discretization of the operator P\f~V 2 P and com- 
pute the relevant matrix element numerically. In doing 
this we can restrict to rotationally-invariant functions, 



meaning we only need to keep track of the radial depen- 
dence. Introducing a radial lattice spacing a we identify 



/1/2 



<H- 



3/2 



5/2 



\ 
I 



f° 



1 



d 
dr 



o 



1 
2^ 



i 1 
-1 1 



V 



(12) 



(13) 



/ 



(note that we begin the lattice half-a-spacing from the 
origin, and that our discretization of d r preserves anti- 
symmetry). The projection operator 



P 



o 



IjVxAT 





(14) 



where the rank of P is related to the size of the hole 
by R — aN. Finally we identify the inner product on 
rotationally-invariant functions 



X 0! (p2 



O 



b\fcj)2 



(15) 



and the state of interest 



d x lx) <-* 




(16) 



(the coefficient is fixed by requiring that the two states 
have the same norm). From these ingredients we can con- 
struct P\ — d r rd r P numerically as an N x N matrix, 



invert it, and take the matrix element between the state 
(16). Extrapolating to the continuum limit N — > oo, 
a — > with R fixed, we find that the charge for a round 
hole is 



Q = 1.28 R 3 



(round hole) 



(17) 



Note that the charge has units of (length) 3 , as could have 
been guessed on dimensional grounds. But one should 
be careful - if two holes merge to form a single hole the 
charges should not be expected to be additive. 



B. Rectangular holes 

In order to gain insight into how the charge depends 
on the shape of the hole, we now consider a rectangular 
hole of size L\ x L 2 . Rather than use a position-space 
lattice, we will work in momentum space with a mode 
cutoff. 

The charge we wish to compute is 

Q= f d 2 xd 2 x' (xKPV^PP)- 1 ^') . 

Inserting complete sets of functions which vanish outside 
the hole, namely 

. . . f ^ T 2^sin(^i) s in(^2) in the hole 

(x\mn) = < vLiL 2 \ hi ) \ l 2 > 

x ' ' [ otherwise 

(18) 
with m, n = 1, 2, 3, . . . we have 



Q 



64LiL 2 



E 



odd m.n,m' .n' 



(mn\(D \m'n') . 

run nn m'lr n'w 

(19) 
Here O = P\f — V 2 P. The matrix elements of O in this 
basis are (assuming odd m and n) 



i\^^^\m'n') 



64Lii 2 



d 2 k 



k\ + k 2 2 cos 2 (k 1 L 1 /2)cos 2 (k 2 L 2 /2)f 1 (m)f2(n}fi(m')f 2 (n') 



(20) 



where 
AM 



k\L\ 



Mn) 



k 2 L 2 



(21) 



In the first line of (20), we used the fact that the projec- 



tion operators do not matter, and in the second line we 
made use of a Fourier integral representation of the basis 
functions (18). 



To compute the charge we introduce mode cutoffs 



m < 2M — 1, n < 2N — 1. Then we can construct 
the operator O numerically as an (MN) x (MN) matrix, 
invert it, and evaluate the sums in ( |19[ ). Extrapolating 
to the continuum limit M, N — > oo then gives the charge. 
For instance, for a square hole with L\ = L 2 = L, this 
procedure gives the charge 

Q = 0.228 L 3 (square hole) (22) 

By comparison, the charge of a round hole ([17]) is 



0.230 (area) 3 / 2 . So a square hole has slightly less charge 
than a round hole of equal area. 

Another simple case is a rectangle with L\ ^> L 2 . In 
this limit the matrix elements ( 20 ) go over to 



i\0\i 



= Sn 




(23) 



i.e., it reduces to a one-dimensional problem. (The easi- 
est way to see this is to rescale k\ — > fci/Zi, k 2 — > k 2 /L 2 
and send L\ — > oo.) Using this in (19) we find that the 
charge for a long thin rectangle is 



Q = In 



/slit 



(thin rectangle) 



(24) 



where Q s iit is the charge for a slit of width L 2 . This is 
something we study in section [TTTt Borrowing the result 
(37), the charge for a long thin rectangle is 



Q = 0.393 L X L\ for! 1 »L 2 



(25) 



Note that, for a given area, a long thin rectangle has less 
charge than a square. This is consistent with the obser- 
vation above, that the value of Q/(area) 3 / 2 is slightly 
smaller for a square hole than for a round hole. It's 
tempting to speculate that round holes have the largest 
charge for a given area. 



C. Square hole: Perturbation theory 



into two parts: a term diagonal in the basis ( 18 1 and a 



nondiagonal term. The nondiagonal part is then treated 
in a perturbation series. 

The fci-integration in p0| ) can be done by moving 
the integration contour slightly above the real k\ axis 
and writing 4cos 2 (fciLi/2) = e iklLl + e~' ifclLl + 2. For 
each term in this decomposition, the integration con- 
tour can be completed in the upper or lower half-plane, 
appropriately. We then pick up pole contributions at 
±(m7r/Li) + ie, ±(m'ir/Li) + ie and cut contributions 
from the square root term \Jk\ + fcf- The pole contribu- 
tions cancel out unless m = ml ' . For the cut contributions 
it is simpler to use the integral representation 



kj 



k 2 



1 



eh. 



H 



" : 2 



fcf + fef + 7 2 



(26) 



and evaluate contributions at the imaginary poles k\ — 

±iy/k 2 +j 2 . 

We follow the same procedure for the A^-integration. 
For a square hole, with L\ = L 2 = L, we then find 



(mn\0\m'ri) = (mn\O po i e + O cnt \m'n') (27) 



where 



(mn\O po i e \m'n') = -\J{m,Tx)' 1 + (nir) 2 5 mm ,S nn > (28) 



The charge Q in ( 19 ) can also be calculated using the 



perturbation method developed in [3] where O is split 



is the diagonal term arising from the real-pole contribu- 
tions in k\ and k 2 integrations, and 



4 
10 



(mn\0,, n \m'n') =-— (rror) 2 / dy vV - 1 (1 + e~™ v ) 



'2^2 



dy I dkk 2 y/y 2 - 1 (1 + er ky ){l + cosk) 
Jo 



(miry) 2 + n 2 7r 2 (rmry) 2 + n' 2 Ti 
nir 



k 2 — n 2 Tt 2 k 2 — n' 2 ir 2 k 2 y 2 + m 2 7r 2 k 2 y 2 + m /2 7r 2 

(29) 



Treating O cut as a perturbation we have 



° - ^polc - pole°cutC pole 



with a corresponding series for the charge, 

Q = Q (0) +Q (1) + ---. 
For a square hole, denoting 

h(m,n) = l/y(r7i7r) 2 + (nir) 2 , 



we find 



(30) q(0) = 64L 3 J2 



odd m.n 



h(m, n) 
(mir) 2 (nir) 2 



0.170 L A 



Q 



(i) 



-64L 3 J2 



(mn\O cn i\m'n'}h(m, n)h(m', n') 



odd m,n.m' ,n r 



(rrnr) (nir) (jn'ir) (n'ir) 



0.039 L A 



(31) 



These first two terms capture, respectively, 75% and 17% 
of the lattice result ( 22 ) , strongly suggesting the validity 



of the perturbation expansion. 















r 





w, 



w- 



A. Lattice approach 



The charge ([33]) can be computed in a variety of ways. 
We begin with a lattice approach, in which we introduce 
a lattice spacing a and discretize the operators via 



dx 2 



(2 1 



<H- 



-1 2 



-1 2 



\ 



(34) 



FIG. 2: Slits of width W\ and W2 and length L, separated by 
a distance r. 



P = 



1-NxN 







(35) 



III. INTERACTION BETWEEN SLITS 

As a simple setting which nicely illustrates our for- 
malism, we conclude by studying a plate with two long 
thin slits in it. We denote the widths of the slits by W 
and the separation between slits by r. We periodically 
identify both Euclidean time with period /3 — > 00 and 
the dimension along the slit with period L — > 00. The 
geometry is shown in Fig. [2| 

The basic result for the interaction energy Q still ap- 
plies, with the following modifications. 

1. The integral over Kaluza-Klein momentum be- 
comes two dimensional, 



BL 



d 2 M 
(2^ 



BL 



[id\i 
~2n 



2. The operator O involves a one-dimensional Lapla- 
cian, and the matrix element (fSp changes to 1 



(xW-ffi x + t P\7>)=~?-K 1 { ii r) 



7rr 



This means that the interaction energy per unit length 
between narrow, widely-separated slits is 



-Bint 

~L~ 



Hdfi ( — #i(jur)) Q1Q2 



1 

An ./,, \nr 

,-3 Q1Q2 



= -5.375 x 10" 



(32) 



where the charge associated with each slit is 

Q= I ' dxdx' (x\(P*/^dlP)~ 1 \x'). (33) 



The rank of P is related to the width of the slit by W 
aN. Finally we identify 



/ 



dx \x) <H- \fa 




(36) 



Computing the matrix element and extrapolating to the 
continuum limit gives the charge for a slit, 



Q = 0.393 W 2 



(37) 



The same coefficient appears in the charge of a long thin 
rectangle (25). 



B. Perturbation theory 

An alternative approach to computing the charge of a 
slit is to use perturbation theory. Inserting a complete 
set of states 



(x\n) 



2 , rmx ^ 

w shl (w 



the charge ( 33 1 becomes 



8W J2 



1,2,3, 



Io->i 



(38) 



(39) 



m, n odd 



where O = P \j —d\P . The matrix elements of this op- 
erator were studied in [5j, where it was found that they 
could be decomposed into 'pole' and 'cut' contributions, 
O = O polc + O cut . The matrix elements for odd m and 



/Tjpolc 

o: 



w 
4 

'nW 



(40) 



cut 

inn 



dyyil 



m 2 ir 2 + y 2 n 2 Tr 2 



1 For d dimensional plates the heat kernel |9l| becomes 
(4ns)~ d ' 2 exp(— r 2 /4s) and the matrix element T8|l becomes 



-2(2^) 2 K(d+l)/2&r). 



2 These are the even-parity modes of [3] equations (29), (32) with 
the dictionary p = m/2, q = n/2, a = W/2. 



The pole contribution corresponds to putting a Dirichlet 
boundary condition on the edges of the slit, while the 
cut contribution captures diffractive effects. Once again, 
treating O cut as a perturbation as in (30 1, the first few 
terms in the series for the charge are 



Q 



(0) 



7C(3) 



W 2 = 0.27UW 2 



Q (1) = 0.0773 W 2 
Q^ = 0.0268 W 2 



(41) 



All terms in this perturbation series are positive. Note 
that the perturbation series appears to be nicely conver- 
gent. Compared to the lattice result (37), the first three 
terms capture 69%, 20% and 7% of t 
spectively, for a total of 96%. 



le exact result, re- 



IV. UNIVERSALITY OF 1/r 7 

The behavior of the potential between two holes in a 
plate is very reminiscent of the van der Waals interaction 
between neutral atoms of zero intrinsic dipole moment. 
In that case, there is an argument due to Feinberg and 
Sucher which shows that the behavior of the potential 
is universal at large distances. In modern language the 
argument for this low energy theorem is that the effective 
action is of the form ^(aEE 2 +aBB 2 ), where E and B are 
the electric and magnetic fields, respectively, and ag, olb 
are the corresponding polarizabilities. If one introduces a 
field $ representing the creation and annihilation of the 
atom as a whole, then the effective action is of the form 



Scff = gid^dvQF^FZ 



92&F 2 



(42) 



For simplicity, we take the atom to be represented by a 
scalar field, the final result is not sensitive to this choice. 



The terms displayed in ( 42 1 are those with the lowest 



dimension; these are the relevant ones for low energy or 
long distance behavior. (The coefficients g\ and gi can be 
related to the polarizabilities, but this is not relevant for 
our argument.) This effective action can be used to calcu- 
late the potential between two atoms at large separation; 



it is given by the two-photon exchange generated by ( 42 1 . 
The computation of this process leads to the 1/r 7 behav- 
ior of the van der Waals potential, in the long-distance 
regime where retardation effects are important [8]. This 
is the universality of the van der Waals interaction. 

A similar low energy theorem applies to the potential 
between two holes on a plate. For this purpose, we first 
note that the boundary action (the action on the holes) 
can be obtained from the quantum effective action T[x] 
which generates the 1PI diagrams. For this, recall that 
the boundary action is defined by 



exp (-S[0o]) - / [#] exp (S[ X + </>]) (43) 

where x obeys the condition x ~ * 4>o as one approaches 
the boundary. We defined % as a special solution of the 



(free) equations of motion with this boundary condition. 
Explicitly, for fields on the left side of the partition of the 
box in which the theory is defined, we could take 



X (x) = / <f-V <j) Q {x') n ■ &G(x\x') 



(44) 



where G(x\x') is the propagator to the left (and right) of 
the plate with Dirichlet boundary conditions. We do not 
need to make this particular choice, \ is an Y specific field 
with the boundary behavior x ~ ** 00 ? so that 4> in (43) 
can be taken to obey Dirichlet boundary conditions. 

Consider now the calculation of r[x] for the left side of 
the partition of the theory in a box. From the Legendre 
transformation of the generating functional for connected 
Green's functions, it is easily seen that we can write 



exp {-T[ X }) = j[d<j>] exp f-S[ X + 0] + J 



ST 
Sx 



<t>) (45) 



where x i s > f° r the moment, an arbitrary field. If we 
choose x to be a solution of the quantum equations of 
motion, ST/5x = 0, obeying the condition x ~ * 4>o ° n 
the boundary, we see that the functional integral on the 
right side of ( [45] ) becomes the defining integral (43) for 
the boundary action. In other words, 



S[4> ]=T[ X ], 



$x 



0, x —* 4>v on the boundary 

(46) 

This gives the boundary action for any theory, includ- 
ing the effect of interactions. It is the quantum effective 
action evaluated on its critical point with the boundary 
condition x ~ * <fio- 

The long-distance interaction between holes comes 
from terms in the effective action which are quadratic in 
X- For these quadratic terms, the low energy or deriva- 
tive expansion of T[x] follows the familiar pattern. For 
a scalar field, the lowest order term is the mass term. 
Since we are considering a massless theory, we can take 
the renormalized mass to be zero, so there is no x 2 term 
in r[x]. The lowest nontrivial term in T[x] is then the 
kinetic term, 



r[x] 



dnxd^x 



(47) 



This term has been canonically normalized by an ap- 
propriate wavefunction renormalization. It leads to the 
boundary action we used above and gives 1/r 7 behavior 
for the hole-hole potential on the plate. Higher dimen- 
sion terms in T[x] will not contribute to the long distance 
behavior of the potential. This argument makes it evi- 
dent that we will get the same behavior for other kinds of 
fields as well. For example, for the electromagnetic field, 
T[A] has similar behavior, with the transverse potentials 
contributing. Hence it will lead to the same 1/r 7 poten- 
tial between holes. It would be interesting to extend this 
result to spinor fields. 

Thus we have obtained a low energy theorem, analo- 
gous to the Feinberg-Sucher proof of universality of the 



van der Waals force: The long distance potential between 
two holes in a plate behaves as 1/r 7 for all massless fields 
for which the leading kinetic operator is the Laplacian. 
(This result holds for a four-dimensional theory.) 



V. DISCUSSION 

In this paper we have analyzed the Casimir interac- 
tion between two holes on a conducting plate (a plate on 
which the fields obey Dirichlet conditions) using the for- 
malism of the non-local field theory developed in [5] . For 
a separation distance r which is large compared to the 
linear dimensions of the holes, the interaction energy is 
proportional to Q1Q2/''' 7 , where Q\, Q2 are charges asso- 
ciated with the holes. The 1/r 7 form is universal, given 
by a low energy theorem. For any number of holes, from 
equation (6j, the result is the same with a pairwise in- 
teraction between holes. (For long slits, because of an 
additional integration along the length of the slits, the 
interaction energy scales as 1/r 6 .) 

There are clearly interesting questions for further anal- 
ysis. The charge Q for each hole depends on the geometry 
of the hole. While we have considered round and rect- 
angular holes and long slits explicitly, it is interesting to 
seek a more general understanding of how Q is related to 
the geometry of the hole. 

The 1/r 7 behavior of the two-hole interaction is for 
large distances. What is the behavior of the interaction 
energy as the holes come closer? While we do not have a 
definitive answer, it is worth noting that we are consider- 
ing the mutual interaction energy between holes; the self- 
energies of the holes do not appear, or, equivalently, they 
have been subtracted out. But when the holes merge, 
the entire energy is the self-energy. It would be inter- 



esting to study the approach and merger of two holes in 
more detail. An intriguing possibility is that, as the holes 
come closer, the potential could reach a minimum value 
at some finite separation, similar to the Lennard- Jones 
potential between neutral atoms. 

A large number of small mobile holes on a plate would 
behave like a gas of particles with a pairwise (two- 
particle) interaction which is attractive and goes like 1/r 7 
at long separations. It would be interesting to study the 
thermodynamics of this "Casimir gas" . 

Our analysis was done using a non-local field theory on 
the holes. There are other methods of analysis, notably 
the world-line formalism, which has been used before for 
studying diffractive effects on Casimir forces [6] . It would 
be interesting, and an independent check, to reproduce 
our results in this formalism. 

Finally, one may ask about experimental observation of 
this interaction. Casimir energies resulting directly from 
fundamental forces, such as the electrodynamic Casimir 
energy, are small and very difficult to measure [9]. How- 
ever the thermodynamic analogue of the Casimir effect, 
which can occur in liquid mixtures near a critical point, 
can be appreciable because, for this effect, temperature 
plays the role of Planck's constant. In fact, there have 
been recent observations of such a thermodynamic ana- 
logue of the Casimir effect 10 . It could be that the 
interaction between holes would be observable in such a 
setting. 
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